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Abstract 

In deriving large n probability distribution function of the rightmost eigenvalue from 
the classical Random Matrix Theory Ensembles, one is faced with que question of finding 
large n asymptotic of certain coupled set of functions. This paper presents some of these 
functions in a new light. 

1 Introduction 

In the study of Edgeworth type expansions for the limiting distribution of the right- 
most eigenvalue from Gaussian Random Matrix Ensembles, we run into finding large 
n expansions of many key functions. For the Tracy- Widom distribution derivation, 
one needs the large n limits of these functions, and they can all be express in terms 
of the couple pair q and p, where q is the Hastings-McLeod solution to Painleve II 
equation behaving at infinity as the Airy function. The frequency of these functions 
in the study of the largest eigenvalue of Gaussian and Laguerre Random Matrix En- 
semble points to the necessity of a study of these functions in their own right. We 
hope this will shed a light into understanding some derivations related to this aspect 
of Random Matrix Theory and related field making use of such functions. 

If one try to read through a proof of an expansion relating various asymptotic 
functions, it's easy to get lost in translation. But if the related functions are well 
known, the reader will probably have a different experience and therefore a better 
understanding of the techniques and tools used for the derivation. 

We present in this paper the derivations of those functions arising in the study 
of the largest eigenvalue for the Gaussian Ensemble of Random Matrix Theory in a 
hope of achieving our goal set above. 

Before stating our results, there is a need to define our functions. 
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For a Gaussian ensemble of n x n matrices, the probability density that the eigen- 
values lie in an infinitesimal intervals about the points Xi < ... < x„ is given 
by 

Fn,p{Xi,- ■ ■ ,Xn) = Cn^exp 1-^ ^x'^A Yl\Xj - Xkf. (1.1) 

V 1 / j<k 

Where /3 = 1 corresponds to the Gaussian Orthogonal Ensemble (G0E„),/3 = 2 
corresponds to the Gaussian Unitary Ensemble (GUE„),and /3 = 4 for the Gaussian 
Symplectic Ensemble (GSE„). 

f ll.ip can also be represented as a determinant involving the variables x^s. For the 
simplest case /3 = 2, we have 

Fn,2{Xl,--- ,Xn) = ^{det[(pj^i{Xi)]f = ^det[Kn,2{Xi,Xj)]ij=i^...^n (1-2) 

nl nl 

with 

V 2 X — y 

and 

1 2/ 

fk{x) = , , Hk{x) ''^ with Hk{x) the Hermite polinomials 

obtained by orthogonalizing the sequence {x'^, k = 0, - ■ ■ ,n — 1} with respect to 
over M. Using this representation, it can be shown that the probability distribu- 
tion function of the largest eigenvalue Xmax is given by the Fredholm determinant of 
the operator with kernel Kn,2 acting on the set {t oo), 

Fn,2{t) = H^max < t) = det(J - i^„,2). (1-4) 

In finding the Edgeworth type expansion of -Fn,2; one needs large n expansion 
of (11.31) or what it amount to, the large n expansion of ipn- In [3], we derived the 
following expression. 

Let the rescaling function r be defined by, 

t{x) = v^2(n + c) + 2"^n"^x, (1.5) 

then 



(/.„(r(x)) = nT. |Ai(X) + i^^-il Ai'(X)n-i + 



[lOc^ - 10c+ ^)XAi(X) 



and 



_ 2 

n 3 



+ 0(n-i)Ai(X)| (1.6) 



^n-Mx)) = |Ai(X) + i^^±ll Ai'(X)n-i + 



(10c2 + 10c+-)XAi(X) 
2 



_ 2 

n 3 



0(n-^)Ai(X)| 



;i-7) 



Ai being the Airy function. These two functions enable us to obtain the following 
expansion of the GUE kernel. 



Kn,2{r{X),r{Y))dT{X) = r' Kn{r{X),r{Y))dX = { Km{X,Y) - c Ai(X) Ai(F)n-3 + 



— {X + Y) Ai'(X) Ai'(F) - + Xr + Y^) ki{X) ki{Y) 
20 

-20c2 + 3 



-(Ai'(X)Ai(F) + Ai(X)Ai'(F)) 



n-i +0(n~^)E(X,r)|rfX (1.8) 

In deriving the finite but large n probability distribution function of the largest eigen- 
value using (11. 8p . and representation (II. 4p we have to factor out of (11.81) the constant 
term (with respect to n) to obtain the representation 



Fn,2{T{t)) = det {{I - KAiiX, r)) . I J + ( J - KAiiX, Y))-' 



cAi{X)Ai{Y)n- 



1 

20 

-20c2 + 3 



{X + Y) Ai'(X) Ai'(F) - (X^ + XY + Y^) Ai(X) Ai(r)+ 



(Ai'(X)Ai(F)+Ai(X) Ai'(F)) 



n-'i +0{n~^)E{X,Y) 



;i.9) 



This Fredholm determinant is computed over the set {t, oo). Thus to complete the 
determination of Fn,2{T{t)) we need to determine the action of the integral operator 
(/ — KpS) on Ai(a;) and Ai(x) where i = 0, 1, ■ ■ ■ . These are the special functions 
in the GUE case, and they are independent of n, they are well known in the literature 
(see for example [lElSllilllSliniEIlEnillHlES]). For these n independent functions, 
we just redefine them here and then introduce their n dependent counterparts. 



K^^X.Y) . A.(A-)Ai(y)-A.(y)A.(AO ^ jy^^^^^ ^^^.^^^^^ ^^^^^ 

p(X, Y- s) = {I- AAi)"'(X, Y; s), R{X, Y; s) = p(X, Y; s) ■ Kp^{X, Y) (1.11) 
this last product is operator multiplication. 



Q,(x;s) = ((/-A^Ai)-\ x^Ai), (1.12) 
Pi{x-s) = ((/-AAi)~\ x^Ai'), (1.13) 
qi{s) = Qi{s;s), pi{s) = Pi(s;s), go(s) := g(s) pq{s) := p{s) (1.14) 



3 



Ui{s) = (Qi,Ai), Vi{s) = {Pi, Ai), Uo{s):=u{s), vo{s) := v{s) (1.15) 

Vi{s) = {Qi,Ai'), Wi(t) = (Pi,Ai'), Wo{s) := w{s), and Vo{s) := v{s). (1.16) 

Here (■,■) denotes the inner product on L^(s, oo) and i = 0, 1,2, These are 
all well known functions, this paper is concerned with the n dependent counterparts 
whose definitions are similar in nature, the changes needed here are on the kernel 
definition. The operator kernel is of the same form as (11.81) 

x-y 



with 



'^(^) = \l^fnix) and ^(x) = ^^ifn-iix) 



Relating this to the previous set of function are the functions ip and ip, they are Ai 
and Ai . We have the following functions. 



pn{x,y;t) := {I - Kn) ^{x,y;t), Rn{x,y;t) := / pn{x, z;t) Kn{z,y;t) dz (1.18) 
these are kernels of integral operators on {t oo) 

POO 

Qn,i{.^]t) ■■= I pn{x,y;t)y'ip{y) dy, Pn,i{x;t) := j pn{x,y;t)y'ip{y) dy (1.19) 



or 

Qn{x; t) := {pn, <f){t oo) Pn{x; t) := (p„, V')(t, oo)- 

And the other functions are 



qnAt) = Qn.,i(^;t), Pn,i(t) = Pn,iit'^'^) ^ qn,o(t) := qn(t) Pn,o{s) := Pn(t) (1.20) 
Un,iit) = iQn,i,^), Vn,i{t) = {Pn,i,^), Unfl{t) := Un{t) , Vnfi{t) := Vn{t) (1.21) 

Vn,i{^) = {Qn,i,i^), Wn,i(t) = {Pn,i,^), Wn,o{t) := w„(t), andVn,o(t) := 5„(t). (1.22) 

Here ( ■ , ■ ) denotes the inner product on L^(t, oo) and i = 0,1,2, ■■■ . 
We will like to point out the following ambiguity in these definitions, the n-independent 
functions have a subscript i whereas the n dependent ones have the subscript n. We 
were not able to find a suitable representations of the set of functions depending on 
the matrix ensemble oi n x n matrices, but the choice of keeping with the original 
Tracy and widom notation was made in part to help the reader go through the topic 
without too much confusion. Thus whenever we use the subscript n we will refer to 
the large size on the underlying matrix ensemble and when i is used it refers to the 
exponent of the variable x appearing in the definition of that specific function and i 
takes values from 0, 1, 2, ■ ■ ■ . One exception is when we will use a second subscript to 
distinguish between the 3 beta ensembles (3 = 1,2,4, in this case we will remind the 
reader of the significance of those values. 
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In deriving the probability distribution function of the largest eigenvalue F„^i(t) 
for the orthogonal ensemble, and FnA{t) for the symplectic ensemble, we encounter 
new sets on functions obeying the same set of relations. 

If we define e to be the integral operator with kernel e{x, y) = |sign(x — y) then 



oo 



Qn,e{^]t) := pn{x,y]t)e{ip){y) dy, qn,e{t) ■= Qn,e{t\t) (1.23) 



oo 



P„,,(x;t) :=y^ Pn{x,y-t)e{ij){y)dy, p„,,(t) := P„,,(t; t). (1.24) 
In a similar way we define 

/oo noo 
Qn,e{x\ t) ip{x) dx, Vn,e{t) ■= J Pn,e{^] t) ip{x) dx. (1.25) 

/OO poo 
Qn,s{x;t) tp{x) dx, and Wn,s{t) = J Pn,e{x;t) tp{x) dx. (1.26) 

And finally we also have for the Gaussian Orthogonal Ensemble 

^n,i(i) := / Rn{x,t;t)dx, Vn,i{t) := / Pn{x;t)dx, Qn,i{t) ■= / Qn{x;t)dx 

J —oo J —oo J —oo 

(1.27) 

(Note here that the second subscript here refers to the beta being 1 for the orthogonal 
ensemble and has nothing to do with the previous discussion on i and n.) For the 
Gaussian Symplectic Ensemble we have, 

/oo poo 
e{x,t)Rn{x,t]t) dx, Vn,4{t) = / e{x — t) Pn{x]t) dx, and 
-oo J — oo 

(1.28) 

/oo 
e{x - t) Qn{x;t) dx, 
■oo 

and the 4 refers to beta being 4 for the Gaussian Symplectic Ensemble. 

We have the large n expansion of most of these functions from previous work. 
What is new in this paper are the large n expansion of Qn,i, Pn,i this can be used 
to derive an expansion for Un,,i, Vn,i, Vn,i, Wn,i- We also have closed formula for 

^n,e; ^n.E) Qn,ey Qn,l, Vn,!, Tln,l, Qn,4, 'Pn,4, and Tln,4- 

In the second section we will give a brief justification of Qn,i and P„ j follow in the 
third section with the justification of these last 9 functions. Again the motivation for 
the derivation of these functions is due to their appearance in the Edgeworth type 
expansion of the largest eigenvalue probability distribution function for the Gaussian 
Orthogonal and Symplectic Ensembles. 
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2 Epsilon independent functions 

Building on O, (OD and we find that 

{I - Kn,2)~\x,y;t) y'ip{y) dy 

/•oo 

:= ((/ - K.,2)-\x, y- 1), y'^{y)) = {I- Ka,2)~\x, y- 1) y'^{y) dy 



therefore we need to find p„(x, t) = {I — i^„,2) ^{x, y; t) in order to find an expres- 
sions for these two functions. But 



(/ - i^„,2)"'(r(X), r(F); r(t)) =/+(/- i^^i)"' (X, F; t) 



cAi(X) Ai(F)r2"3- 



1 

20 



(X + Y) Ai'(x) Ai'(r) - (x^ + xy + r^) Ai(x) Ai(y)+ 



-20c^ + 3 



(Ai'(X)Ai(F)+Ai(X)Ai'(y)) 



/ + cQ{X) Ai(F)n-3 - - 



{Q^iX) + Fgi(X) + F2g(X)) Ai(F) 



+0{-)E{X,Y) 
n 



n-r^ +0{n-^)E{X,Y) 

{Pi{X) + YP{X))M{Yy 
20c2 + 3 



I - Km{X,Y))-' 
(2.1) 



-(p(x)Ai(r) + g(x) Ai'(F)) 



n 3 



/ - cQ{X) Ai(F)n-3 + 



20 



I - Km{X,Y))-^ 



(Pi(X) + rP(X))Ai'(F)- 



(Q2(X) + YQ,{X) + F2g(X)) Ai(F) + ^ (P(X) Ai(F) + Q(X) Ai'(F)) 



-20c2Q(X;s)m(s) Ai(y) 



1 



n-3 +0(-)E(X,F) 



n 



|-(/ - J^Ai(X,F))-^ = p„(X,F;t) 



Note that with this representation of p all the Qn,i and P„^j will have no term inde- 
pendent of n, but only Q„ o •= Qn and P„ q := Pn, in [3] we find that 



Q„(r(X);r(s)) =n6 



g(X;s) 



2c- 1 



P{X-s)-cQ{X-s)u{s) 



n 3 



(lOc^ - 10c + ^)Qi{X- s) + P2(X; s) + (-30c2 + 10c + ^)Q{X- s)v{s) 
-Pi(X; s)v{s) + P(X; - QalX; s)n(s) - gi(X; s)ni(s) - Q(X; 5)^2(3) 



-lOc^ + ^)P{X; s)u{s) + 20c^Q{X; s)u\s) 



n 3 



20 



0{n-^)E,{X-s) 



(2.2) 



and 



P„(r(X);r(s)) =na 



Q{X-s) + 



2c + 1 



P{X-s)-cQ{X-s)u{s) 



n 3 



+ 



(10c' + 10c + ^)Qi{X- s) + P2{X- s) + (-SOc' - 10c + ^)Q{X- s)v{s) 



-Pi(X; s)v{s) + P(X; s)v^{s) - Q^iX- s)u{s) - Qi{X; s)ui{s) - Q(X; 3)^2(5) 



-lOc' + ^)P{X; s)u{s) + 20c^Q{X; s)u'^{s) 



_ 2 
n 3 



Oin~')Ep{X-s) 



(2.3) 



Using 



<3„,(t(A').t(3)) = (p„(T(A'),T(y),T(s)),(T(r))V(T(F)))M,) „o) 



E 

/c=0 



i\ 22 (n + c) 2 



(p„(r(X), r(F); r(.)), FV(r(r)))(,(.) 00). 



and 



XV(r(X)) = riH <j X^Ai(X) + i^^^X^ Ai'(X)n' 



(10c2-10c+^)X^+iAi(X) 



+ X^+2Ai'(x) 



_ 2 

n 3 



+ 0(n-^)Ai(X) 



(2.4) 



we find that 



p(X, Y- s)-^{t{X)) = nt { QkiX; s) + ^^^-^P,(X)n-i + 



(10c2-10c+-)Qfc+i(X) 



+ Pk+2{X) 



n 3 



20 



0{n-')QkiX) \. 



(2.5) 



Combining this with the action of the first factor on the right of (12. ip gives the 
following expression for (p.„(r(X), r(F); r(s)), X'^(y9(r(X)) 



n-e \QkiX;s) + 



2c- 1 



Pk{X-s)-cUk{s)Q{X-s) 



n 3 



^ + [(10c-20c2)5fc(s)g(X;s) + 



(10c'-10c+-)Qfc+i(X;s)+Pfc+2(X;s)+Pi(X;s)i;fc(s) + P(X)(rAi (y),Qfc(F))(, 00) 
-Uu{s)Q2{X- s)-Q,{X- s){Y Ai(F), Qfc(F; s))(, oo)-Q(^; s)(>^' Ai(F), gfc(r; s))(, 00) 



—OOr"^ _l_ 3 — 20c^ -I- S 
+ TT^PiX; s)ukis) + ^Q(X; s)v,is) + 20c2g(X; s)w(s)w,(s) 



2 

n 3 



-Oi-)EiX;s)}. 



n 

To simplify the inner product in this last expression, we use the following recurrence 
relation derived in [21] 

Qk{X; s) = X^Q{X; s) - T.i+j=k-i;i,j>oi'"jQi - to have 

noo 
X Ai(X)p(X, Y; s)Y' Ai(F) dV dX 

= (Qi(X;s),X^Ai(X)) = {XQ{X;s) + u{s)P{X;s)-v{s)Q{X;s),X''Ai{X)) = 

Uk+i{s) + u{s)vk{s) - v{s)ukis) 

and 

(Qfc(X;s),X2Ai(X))(, ^) = (Q2(X;s),X^Ai(X))(, ^) = 
= {X^Q{X- s) - v{s){XQ{X- s) + u{s)P{X; s) - v{s)Q{X- s)) - u{s){XP{X- s) 
-w{s)Q{X- s) + v{s)P{X- s)) - v,{s)Q{X; s) + ni(s)P(X; s),X'' Ai(X)) 
= Uk+2{s) - v{s)uk+i - v{s)u{s)vk{s) + visfukis) - u{s)vk+iis) 
+u{s)w{s)uk{s) - u{s)v{s)vk{s) - vi{s)uk{s) + ui{s)vk{s). 

we also have 

(Q,(X;s),XAi'(X))(, = {Pi{X;s),X''Ai{X)) = 
Vk+iis) + v{s)vkis) - w{s)ukis). 

We therefore have 



g„,(r(X),r(s)) = ^ 
2c- 1 



i\ 22 ''(n + c)~ 



k=0 



k\{t-k)\ n 



2 6 



.\Qk{X-s)+ 



-Pk{X-s)-cuk{s)Q{X-s) 



n 3 + 



(10c-20c')Cfc(s)g(X;s) + 



(10c2-10c+^)gfe+i(X; s)+Pk+2{X- s)+Pi{X- s)vk{s)+P{X) {vk+i{s)+v{s)vk{s)-w{s)uk{s] 
-Uk{s)Q2iX; s) - Qi{X; s) I Uk+i{s) + u{s)vk{s) - v{s)uk{s] 



-Q{X; s) Uk+2{s) - v{s)uk+i - v{s)u{s)vk{s) + v{s) Uk{s) - u{s)vk+i{s) 



+u{s)w{s)uk{s) - u{s)v{s)vk{s) - vi{s)uk{s) + ui{s)vk{s) 



—OOr'^ -I- 'X — 9nr^ -I- '\ 

+ , P(X; s)uk{s) + g(X; + 20c2g(X; 



_ 2 

n 3 
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+0{-)E{X;s)\. 
n J 

In a similar way we have 

P.,(r(X),r(.)) = (p„(r(X),r(F),r(s)),(r(r))V(r(F)))M,) 
i\ 25(n + c)^ 



(2.6) 



fc=0 



k\ (z - A;)! 2^n 



-(p„(r(X), r(F); r(s)), F'^V^(r(F)))(.(,) 



And (p„(r(X),r(F);r(s)),X'=V^(r(X)) is equal to 
■2c + 1 



Q,(X;s) + 



n 3 



5 + [-(10c + 2Q^)vk{s)Q{X- s) + 



(10c2 + 10c+-)gfc+i(X; s) + Pfc+2(^; s) + Pi(X; sK(s) + P(X)(F Ai (F), Qfc(F))(, oo) 
-Uu{s)Q2{X- s)-Qi{X; s){Y Ai{Y), Qk{Y; oo)-Q{X; s){Y' Ai{Y), Qk{Y; oo) 

-20c2 + 3 



-20c2 + 3^,^ , , , 
+ ^ P{X;s)ukis) 



-Q{X; s)vk{s) + 20c"Q(X; s)u{s)uk{s) 



_ 2 

n 3 



20 



-Oi-)E{X;s) 
n 



this therefore gives 

i 

P„,(r(X),r(s)) = ^ 



22 '^{n + c) 2 



fc=0 



k\ {i - k)\ ni 



k_l 
6 



Qk{X-s) + 



2c + 1 



Pfc(X;s)-cnfc(s)Q(X;s) 



n 



+ [-(10c + 20c')vk{s)QiX; s) + 



(lOc^ + 10c + -)Qfc+i(X; s) + Pfc+2(X; s) + Pi(X; s)i;fc(s) + 



P(X) yVk+l{s) + v{s)Vk{s) - w(s)Mfc(s; 



Uk{s)Q2{X; s) - Qi(X; s) ( Uk+i{s) + u{s)vk{s) - v{s)uk{s) 



-Q{X; s) (^Uk+2{s) - v{s)uk+i - v{s)u{s)vk{s) + visfukis) - u{s)vk+i{s) 



+u{s)w{s)uk{s) - u{s)v{s)vk{s) - vi{s)ukis) + ui{s)vkis] 



_20c^ -I- 3 — 20c^ -I- 3 
+ TT^PiX; s)uk{s) + ^^Q(X; s)vkis) + 20c^Q{X; s)uis)u,is) 



_ 2 

n 3 
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-Oi-)EiX;s)\. 



n 



(2.7) 



When we set i to zero we recover s) and Pn{X\ s). 

We see immediately that these two series representations of Qn.i and P^^i are not 
in terms of when i is not zero. We can use (12.51) . (II. 7p . (12. 6p and (12. 9p . to derive 
an expansion for Un,i, Vn,i, Vn,i and Wn,i from their representations 

UnA't) = (Qn.,i(a;;t),V5(x))(t oo) VnAt) = {Pn,i{^\t) , ^{x)) {t oo) 

Vn,i{^) = ((5n,i(a:;t),^/'(x))(t oo) and Wn,i{t) = {PnA{x;t),tp{x))(^t oo)- 



We would like to note that (12. 6p and (12.90 are the new quantities in this section, as 
additional corollary the derivation of qn,i{t) and Pn,i(t). 



9n,i(r(s 



)) = E 



k=0 



k\it-k)\ n 



k_l 
2 6 



• i qkis)+ 



2c- 1 



Pk{s) - cufc(s)g(s) 



n 3 + 



{10c - 20c')vk{s)q{s) + 



(lOc^ - 10c+ -)gfc+i(s) +pfc+2(s) +pi{s)vk{s) +p{s) \^k+iis) + v{s)vk{s) - w{s)ukis 

-Uk{s)q2{s) - qi{s) (^k+i{s) + u{s)vk{s) - v{s)uk{s) 

-q{s) (^Uk+2{s) - v{s)uk+i - v{s)u{s)vk{s) + v{sYuk{s) - u{s)vk+i{s) 

+u{s)w{s)uk{s) - u{s)v{s)vk{s) - vi{s)uk{s) + ui{s)vk{s] 
-20c2 + 3 



H — ^'^'l '^^ p{s)ukis 



-Q{X; s)vk{s) + 20c q{s)u{s)uk{s) 



n 3 



20 



+0(-)e,(.) 



77, 



(2. 



Pn,i(r(s 



)) = E 



22"''(n + c)~ 



k=0 



2 6 



■X 



gfc(s) + 



2c + 1 



Pk{s) - cuk{s)q{s) 



n 3 



^ + [-(10c + 20c2)0fc(s)g(s) + 



(lOc^ + 10c + -)qk+i{s) +Pk+2{s) +Pi{s)vk{s) 



p{s) Vk+i{s) + v{s)vk{s) - w{s)uk{s] 
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Ukis)q2is) - qi{s) + u{s)vkis) - v{s)uk{s] 



-q{s) Uk+2{s) - v{s)uk+i - v{s)u{s)vk{s) + - u{s)vk+i{s) 



+u{s)w{s)uk{s) - u{s)v{s)vk{s) - vi{s)uk{s) + ui{s)vkis] 



_ 2 

n 3 

~2(r 



+0{-)ep{s] 



n 



(2.9) 



In [5], we found an expression for Rn{x,y) = Pn{x,y) ■ Kn^2{x,y), this also follows 
from fll.Sp and fl2.ip . Note that the following representation will give the same result, 



Rnix,y;t) 



Qn{x; t)Pn{y; t) - P^ix; t)Qn{y; t) 



x-y 

RniriX), r(r); r(s))rfx = R{X, Y; s) - cQ{X; s)QiY; s) n 

_ 2 

n 3 



(2.10) 



+ 
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Pi(X;s)P(r;s) + P(X;s)Pi(r;s) 



-Q2{X- s)Q{Y- s) - Qi{X- s)Qi(Y; s) - Q{X; s)Q2iY; s) + 20c\o{s)QiX- s)Q(Y- s) 
3 - 20c2 



{P{X;s)Q{Y; s) + Q{X; s)P{Y; s)) 



0{n~')en{X,Y) 



dX. (2.11) 



3 Epsilon dependent functions 

The corresponding epsilon functions come from the study of the leftmost eigenvalue 
from GOE and GSE. We present here the system of equation satisfied by those func- 
tions and a solution to these equations leading to our desired functions. 

To simplify notations we define 

Vn,e{t) = 1-Vn,e{t), and (t) = 1 - 7^„,l (t) . (3.1) 

With this notation, system is 

VnAt) = Pn{t) ■ Vn,e{t) \ , (3.2) 

qn,e{t) J \ -Pn{t) g„(t) / V g„,,(t) 

the boundary conditions in this case are 

K,.(oo) =1 1 I . (3.3) 
qn,e{oo) J 
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For the orthogonal ensemble 



Qn,l(t) \ / qnit)\ f Qn,l{t) 

'Pn,l{t) = Pn{t) ■ Vn,l{t) | , (3.4) 



with boundary conditions in this case are 



Qn,l(00) ^| 




f 2c, 


■Pn,l(00) 







^n,l(00) J 







as n is even. (3-5) 



We also have for the symplectic ensemble 



VnAit) = -Pn{t) ■ VnAit) I . (3.6) 

nnAt) ) \ -Pnit) -qnit) / \ n^^^t) 

where IZnAif) = 1 + 'Tln,Aif), with corresponding boundary conditions 

-c^ \ / 




-Ctp = I — c^) I as 77, is odd. (3.7) 

The first two set of equations were solved in [3] , here we give the general solution from 
the series expansion derived there. We will not go back into the derivation, but would 
like to point out that this is the direct consequence of those matrix exponentials. Our 
goal here is to give a close formula for those functions. 
We define 

Qnix) dx and b{t) = / p„(x) dx. (3.8) 

We note that these two functions scale (under the transformation r) in the large n 
limit to the same function 

1 f°° 1 

We give this to justify our notation used bellow, and it says that for very large n, the 
argument of all the hyperbolic functions is real. With this notation, we have 



UnAt) = - VMt)Kt)] + ^ si^h VMtW), (3.9) 

VnAt) = + cosh VMMt)] - c^J ^ sinh V2a(t)6(t), (3.10) 



or 



Vr, 



^{t) = 1 - ^[1 + cosh ^2a{t)b{t)] + PjK sinh ^2a{t)b{t) (3.11) 
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and 



qn,e{t) = -J ^ sinh ^2a{t)b{t) + cosh ^2a{t)b{t). (3.12) 



This result is vahc^ for the GOE when ^ and for the GSE we have = 
In the same way we find that for the GOE case, the calhgraphic functions are 



Qn,iit) = c^[l + cosh ^2a{t)b{t)] - J ^ sinh ^2a{t)b{t), (3.13) 



Vn,i{t) = c^^[cosh ^2a{t)b{t) sinh ^Mfm, (3.14) 



and 



or 



nn,i{t) = -2c^ J^ly sinh ^2a{t)b{t) + cosh ^2a{t)b{t) (3.15) 



nn,i{t) = 1 + 2c^J^^sinh ^/2a{t)b{t) - cosh ^2a{t)b{t) (3.16) 



where 



c, 



.^-y7:nf"2-^"-'''i^—. (3.17) 

nj Z). 



A large n expansion for w„ g„ ^ is given in |1] on page 17. A large n expansion of Vn^i 
is equation (3.58) and 7ln,i is equation (3.59) of the same work. We will therefore give 
here an expression for Un,e and Qn,i for large n. Substitution of a{t) = qn{x) dx 
and b{t) = pn{x) dx into (13.91) and f l3.13p yields the following results. 

Theorem 3.1. for s bounded away from minus infinity, 

UnAr{s)) = ^(l-e"^(^)) + (^^(e~^(^) + cosh(/i(s))-2)-^e-'^(^)^n-5 + 

(e-''(^) {uisf (- (-1 + e^(^)) (-5 - 12c + (3 + 4c)e'^(^)) - 
2/i(s) (1 + 6c - 2ce2'^« + 4cV(s))) + 

q[x]u[x] dx+ 

/ / POO PCX) \ 

8/i(s) -lOc (-3 + e'^^")) (-1 + e''^")) / [a;] dx — qi[x] dx ) + 



fi{s) ({3 -20c') £ 



p[x]n[a;] dx + 3 q[x]v[x] dx+ 



^The computation for GOE assumes n to be even and for GSE assumes n to be odd 
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f [x]pi[a;] + 2 J p2[x]dx + 3 J qi[x]dx — c y{j q[x]u[x]dx) — 

30 POO POO \ \ / POO 

20 ( 2 / dx — 3 q[x]v[x] dx + qi[x]dx] ] — 2 i u[x]q2[x] dx+ 

qi[x]ui[x] dx + / q[x]u2[x]dx— / c/x ) ) ) ] n^^ + 0{n^^) 



We also have 



VnAris)) = J(l-e-'^W) + (^sinh(/i(^)) + ^e-'^(^))n4 



'Ms) 

1 

+i^(s)^(cosh[/i(s)]/i(s) (-1 + 4c - 4cV(s)) + (l - 4c + + 4cV(s)^) sinh[/i(s)]) 



4cz/(s) / dx (- cosh[/i(s)]/i(s) + 2ce ''^"^ fi{s)'^ + sinh[/i(s)]) 



4/i(s) |^e~''(^V('5) (^(-3 + 20c^) y fix- 3 y g[x]i;[x] rfx-2 J v[x]pi[x]dx 

/OO /'OO / POO / POO 

P2[x] dx — 3 qi[x]dx + i { q[x]u[x] dx^ — 20 I 2 / g[x]n[x]^ dx 



OO POO \ \ / POO POO 

3 I q[x]v[x\dx + / gi[x](ix) 1 + 2 1 / 'u[x]g2Nc?x+ / gi[x]Mi[x] c/x + 



g[x]M2[x] dx — I p[x]t'i[x] dx]j + 20c g[x]t'[x] dx — j qi[x\ dx ) sinh[/i(s)] ) \ n 



we also have 



Iff f°° 

4cz/(s)( / g[x]M[x] dx) (/i(s) (cosh[/i(s)] + 2ce"^^^V('S)) - smh[/i(s)]) + 



8V2fi{sy 

iy{s)^{cosh[i2{s)]i2{s) (l + 4c + 4c^;u(s)) - (l + 4c + + 4cV(s)^) sinh[/i(s)]) + 

/ / poo poo poo 

4:fi{s) (e~^^''^n{s) i {-3 + 20c^) / ^»[x]n[x] rfx- 3 / g[x]w[x] dx-2 / t;[x]pi[x] dx 



J s 

OO fOO / POO / POO 

2 I f I „r^l„.r^l on I n I „r^l„.r™l2 



2 I p2[x]dx — 3 J qi[x]dx + c g[x]n[x] dx) —20^2 J g[x]n[x] dx- 

POO POO \ \ / POO POO 

3 I q[x]v[x]dx+ I qi[x]dx \ | + 2 1 / ti[x]g2[2^] dx + / gi[x]ui[x] dx 

p[x]fi[x] dx 
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20c y— J q[x]v[x]dx + J qi[x] dxj sinh.[jj,{s)]j j n 3 

and the GOE„ calligraphic variables are 
Theorem 3.2. for s bounded away from minus infinity, 



-4ci^(s) f q[x]u[x] dx (cosh[yu(s)] + 2ce~'^^''^ n{s)) - 



8\/2/x(s)- 

sinh /i(s))+i^(s)^(cosh[/i(s)]/i(s) (l + 4c + 4cV(s)) - (l + 4c + ^(s) + 4cV(s)^) sinh fi{s)) 



4/i(s) |^e"''^*V('5) (^(-3 + 20c^) y j9[x]n[x] dx- 3 y q[x]v[x]dx-2 J v[x]pi[x]dx 

/•OO /"OO / /»oo / /»oo 

—2 / p2[x]dx — 3 / c/x + c^ I ( / (ix)^ — 20 I 2 / g[x]M[x]^(ix 



00 



-3 J q[x]v[x]dx + J qi[x] dx j j + 2 u[x]q2[x]dx + J qi[x]ui[x] dx+ 

o /-OO \ \ / POO roo \ 

g[x]u2[x] dx — p[x]vi[x] dxj \ + 20c I / qi[x] dx — q[x]v[x] dx j sinh/i(s) ) ) 
we also have 

((e~^('^ {u{sf ((-1 + e^W) (5 - 12c + (-3 + 4c)e^(')) - 



16 {V2ii{sf) 

2(/i(s)) (1 + 2c (e^^^^) - 3) + 4cV(s)))+4cz/(s) (4e^(^) - 3 + e^^'^'\^l{s) - 1) - 3/i(s) + Ac^l{sf) 

/ POD POD \ 

q[x]u[x] dx+8fi{s) (lOc (-3 + e^^'^) (-1 + e^^'^) ( / q[x]v[x] dx - qi[x] dx ) + 



I (3 — 20c^) / p[x]m[x] c/x + 3 / g[x]f [x] (ix+2 / t>[x]pi[x] (ix+2 / p2[x]dx+ 

\ J s J s J s J s 

/ POO / POO POD 

qi[x] dx — i { g[x]n[x] dx)"^ — 20 I 2 / g[x]n[x]^ dx — 3 g[x]f [x] dx+ 



OO 



2 



qi [x] dx 

q[x]u2[x]dx— I p[x]vi[x]dx] ) ) ) ) n~3 
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and the last of the GOEn function is 

7e„,i(r(s)) = (l-e-'^W) + 
1 



iy{s 



Ms. 



sinh fi{s) — cq{s)e 



n 3 + 



^/i(s) 



{{Acu{s)fi{s) (- cosh[/i(s)]/i(s) + 2ce~'''^'^fi{sf+ smh[/i(s)]) + 



i/(s)^(cosh[/i(s)]/i(s) (-1 + 4c - 4cV(s)) + (l - 4c + /i(s) + 4cV(s)^) sinh[/i(s)]) - 

/ / POO POO POO 

4/i(s) ( e"^('V(-5) ( (-3 + 20c^) / p[x]u[x]dx- 3 q[x]v[x]dx-2 v[x]pi[x]dx 



-2 / p2[x]dx — 3 / c/x + c^ I ( / g[ 

Js J s \ J s 



dxf - 20 2 / 

u[x]q2[x\dx+ / gi[a;]ni[a;] 
g[a:]u2[x] dx — I p[x\vi[x\ dx \ \ + 

20c q[x]v[x]dx — J qi[x] dx^ smh[fi{s)]^^ n~ ^ + O {n~^) 

For the symplectic ensemble, the calhgraphic variables were mentioned in the 
similarity of the corresponding system of differential equations to the GOE calli- 
graphic system makes our derivation simpler, the coefficient matrices are the same 
up the minus sign. Keeping with the same notation, we see that 



llnAit) 



/ 1(1 + cosh V2^) f (cosh V2a6- 1) ^sinhV2a6\ 
^(coshA/2a6- 1) 1(1 + cosh A/2a6) a /j^ sinh A/2a6 







V V ^ ^vah y2ab ^J~^ sinh y/2ab cosh \/2ab y 



(3.18) 

We dropped the t dependence of a and b in the above matrix for esthetic reason. This 
gives 



Qn,4(t) = -C^^[cOsh y^MtWt) - 1] 



a{t) 



2b{t) 



sinh v/Mt)Kt) , (3.19) 



K,4(t) = -c^^[l + cosh ^2a{t)b{t)] - J ^ sinh .J2a{t)b{t), (3.20) 
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and 



or 



n^At) = -c^J ^ sinh ^/2a{t)b{t) + cosh ^2a{t)b{t) (3.21) 



TZnAt) = -c^,y ^ sinh y/2a{t)b{t) + cosh ^/2a{t)bit) - 1. (3.22) 

for the GSE, we have the corresponding formula for m„ e, Vn,e and qn^e- 
Theorem 3.3. for s bounded away from minus infinity, 



Un,e{r{s)) = -sinh (^-) + 1^^^;]^'^^°^^'^/^'^^)) ~ ^) ^ sinh(/i(s)) ) n 3 + 

/■oo 

8cz/(s) / (ia; (— 1 + cosh[/i(s)] (l + c/i(s)^) - 



16/i(s)2 

/i(s) sinh[;u(s)])+i^(s)^(4+8c-4 cosh[/i(s)] (l + 2c + cV(-5)^) + (1 + 8c)yu(s) sinh[/i(s)]) + 



4/i(s) ^40c(— 1 + cosh[/i(s)]) ^— y [x] c/x + y + 

/ POO / POO 

/i(s) I — cosh[yu(s)]( / q[x]u[x]dxy+ I (— 3 + 20c^) / p[a;]n[x] rfx— 

1 /"OO /»00 POO 

q[x]v[x]dx — 2 / t>[x]pi[x] (ix — 2 / p2[x]dx — 3 / gi[x](ix— 

J s 

) ;"cxi /"oo \ / POO 

20c^ ( 2 / q[x]u[x]'^ dx — 3 / g[x]f [x] dx + / gi[x](ixJ+2( / n[x]g2[2;] (^3;+ 



gi[x]ni[x] dx + / g[x]u2[x](ix 

+0(n-i) 



J p[x]fi[x] (ix^^ sinh[/i(s)]^^^ n s 



and 



2//^(s)x , cq{s) 



Wn,e(i"(s)) = -sinh (^— ) H ^sinh(/x(s))n 3 + 



1 



4c cosh[/i(s)]/i(s)(i/(s) — / g[x]M[x] (ix) + 



16/i(s) 

z/(s)^ + 4/i(s) ( (-3 + 20c^) / p[x]u[x]dx-3 / g[x]t'[x]rfx- 



oo 



2 / [x]pi [x] dx — 2 J p2[x]dx — 3 J qi[x]dx — 20c y2 J g[x]M[x] dx- 

30 POO \ / POO POO 

q[x]v[x]dx+ / gi[x](ix)+2( / M[x]g2N c^a^ + / gi[x]Mi[x] (ix+ 
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— J p[x]vi[x] dx^^^ smh.[fi{s)]^ n 3 + 0(n ^] 



we also have 



1 f i^(s) . , , , cq(s) 



gn,£(r(s)) = --^sinh(/i(s)) + (^^-^=--^ sinh(/i(s)) + — ^ cosh(/i(s)) ) n 3 + 



1 



cosh[/i(s)]/i(s) I (1 + 4c)i^(s) — 4cz/(s) / g[x]n[a;] c/x+ 



4yu(s) I (— 3 + 20c^) / p[x]u[x]dx — 3 / q[x]v[x]dx — 2 / v[x]pi[x] dx — 2 



P2[x]dx—3 / (ix— 20c^ I 2 / g[x]-u[x]^ (ix — 3 / (ix + / gi[x](ix) + 

Js V Js J s J s 

/*00 1*00 l*QO 1*00 \ 

2( / M[x]g2[a^] c^a^ + / gi[x]Mi[x] (ix + / q[x\u2[x]dx — I p[x\vi[x]dx 

POO 

(1 + 4c + 4cV(s)^) -4cz/(s) (l + 2c/i(s)^) / q[x]u[x]dx + 

J s 

POO / POD POD \\\ 2 

4c/i(s) ( c/i(s)( / g[x]n[x] dx)^ + 20 I / g[x]t4x] rfx — / gi[x](ix) ) ) sinh[yu(s)])?7,~3 



+0(n-i). 

For the GSE Calligraphic functions we have the following expansions, 
Theorem 3.4. for s bounded away from minus infinity, 

Qn,4(T(s)) = -^-= (1 - cosh/i(s) + 2 sinh /i(s)) + 
2v 2 



2v^/i(s) 2v^ 
1 



1) ^(2 cosh/i(s) — sinhyu(s)) J n 3 + 

2v2 / 

(e-^(^) {u{sY (-2 (-1 + e^(^)) (-3 - 8c + e^^^^) - 



32V2/i(^ 

(3 + 16c + e2^(^)) + 4c2 (-3 + e^^^^^)) fi{sf) - 8cu{s) 

/OO 
q[x]u[x] dx + 

(/ POO POD \ 

80c (-1 + e'^^")) I / g[x]t4x] dx - qi[x] dx J + 

POO POO 

(-3 + 20c2) (3 + e''^^^)) / p[x]n[x] dx + (-3 + e^^^^^) ( / g[x]n[x] dx)''+ 

J s J s 

(POO POO 
AOc^ / q[x]u[x]^ dx+{3- 60c^) / g[x]f [x] dx+ 



2 / da; + 2 / p2[x] dx + [3 + 20c^^ / qi[x]dx — 2i / n[x]g2[a;] 



oo 



gi[x]Mi[a;] (ix + J q[x]u2[x] dx — J p[x]vi[x] dx j j j j j n 3 + 0(n ^) 
the next function is 

K,4(7-(s)) = ^^(2sinh/i(s) - cosh /i(s) - 1) 

+ (^ 2^^f!{s) ^^^^ " ^^^^ ~ ^^^^ ^'^^''0 '^'^ 

1 /""^ 

((8cz/(s) / cosh[/i(s)]/i(s) + 



16V2/i( 

c/i(s)^(cosh[/i(s)] — 2 smh[/i(s)]) + sinh[/i(s)]) 
i/(s)2 (-2 cosh[/i(s)]/i(s) (1 - 4c + 2cV(s)) - 
(2 - 8c + + 8cV(s)^) sinh[/i(s)]) - 



"OO 

4/i(s) I yu(s) I c^( / (ix)^(cosh[/i(s)] — 2 sinh[/i(s)]) + 

(-3 + 20c2) ( / p[x]u[x] da;)(2cosh[/i(s)] - sinh[^(s)]) - 

J s 

roo 1*00 /»oo 

40c^ / q[x\u[xY dx + {3 — QQc^) I q[x]v[x]dx + 2 j v[x\pi[x]dx+ 



u[x\q2[x\dx+ I qi[x\ui[x\ dx- 



1*00 nOQ / nc 

2 / p2[x\dx+ [?> + 2Qc^) I qi[x\dx-2{ j 

J s J s \J s 

q[x]u2[x] dx — J p[x]vi[x] dx^^ (2 cosh[/i(s)] — sinh[/i(s)])^ + 
J qi[x] dx^ sinh[/i(s)]^ ^ n~3 + 0{n~^) 



40c yj q[x]v[x] dx 
the last of our function is 



7?,„^4(r(s)) = ^cosh/i(s) — - smh/i(s) — 1 ) + 

vis] CQ is] 

sinh H — (cosh — 2 sinh fi{s)) ) ri~3 



4/i(s) ' 2 

1 



((— 4ci/(s)( / (ix)(cosh[/i(s)]/i(s) 



16/i(s)2 

(1 + 4c/i(s)) - (1 + 2c/i(s)2) sinh[/i(s)]) + 
i/(s)2 (cosh[/i(s)]/i(s) (1 + 4c + 8cV(s)) - 
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(l + 4c + 2/i(s) + 4cV(s)^) smh[/i(s)]) + 4/i(s) 
|^(-3 + 20c^) p[x\u[x\ rfa;(cosh[^(s)] - 2 sinh[^(s)])- 

/•OO /"OO /"OO 

40c^ / (ix + (3 — 60c^) / [x] (ix + 2 / t>[x]pi[x] 



H c^a; + (3 + 20c^) y qi[x]dx — 2yj u[x]q2[x]dx + J qi[x]ui[x] dx+ 
q[x]u2[x]dx— / dx 1 ) (cosh[yu(s)] — 2 sinh[yu(s)]) + 



c ( J q[x]u[x]dx) (2 cosh[/i(s)] — sinh[/i(s)])J + 
20c ( — / g[x]f [x] dx + qi[x] dx ) sinh[/i(s)] ) ) n~3 + 0{n^^) 
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